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MTI  NOISE  INTEGRATION  LOSS 


INTRODUCTION 

MTI  signal  processing  correlates  the  receiver  noise  and  thus  results  in  degraded  detec- 
tion performance  when  the  MTI  pulses  are  integrated.  Previous  investigators  [1,2]  have 
described  the  decreased  performance  in  terms  of  a reduction  in  the  effective  number  of 
independent  pulses  integrated.  However,  since  the  effective  number  of  pulses  can  be 
represented  by 


N„  = . 

(o2/m2)„y^ 

where  a and  m are  the  standard  deviation  and  mean  of  the  input  samples,  Ng  has  a precise 
meaning  (in  terms  of  detection  performance)  only  if  the  output  noise  distribution  is  com- 
pletely specified  by  Ng.  For  instance,  when  the  number  of  pulses  integrated  (N)  is  large, 
the  integrated  output  is  approximately  Gaussian  distributed  and  integration  improvement 
varies  as  the  square  root  of  the  number  of  pulses  integrated.  Thus  the  loss  (due  to  the 
MTI  correlating  the  receiver  noise)  in  signal-to-noise  ratio  (S/N)  for  a large  number  of 
integrated  pulses  is 

L = 10  log  (N/Ng)^'"^. 

In  this  report  the  MTI  integration  loss  is  calculated  when  the  number  of  integrated 
pulses  is  small  and  thus  the  output  is  not  Gaussian  distributed.  This  calculation  is  per- 
formed using  simulation  techniques.  First,  the  appropriate  thresholds  for  a given  probabil- 
ity of  false  alarm  Pfg  are  calculated  using  importance-sampling  techniques.  Next,  proba- 
bility of  detection  Pp  curves  are  generated  by  simulation  of  the  pulse-by-pulse  video. 
Finally,  the  MTI  integration  loss  is  found  by  comparing  the  generated  Pp  curves  with  those 
for  independent  samples  [3]. 


FALSE-ALARM  THRESHOLDS 

Although  Monte-Carlo  simulations  have  been  used  for  many  years  to  calculate  Pp 
curves,  they  have  not  been  used  to  calculate  Pfg  curves  because  of  the  enormous  numl^er 
of  repetitions  usually  required;  approximately  lOjPfg.  However  this  difficulty  can  be 
overcome  by  using  importance  sampling  [4].  The  fundamental  principle  of  the  importance- 
sampling technique  is  to  modify  the  probabilities  that  govern  the  outcome  of  the  basic 
experiment  of  the  simulation  in  such  a way  that  the  event  of  interest  (the  false  alarm) 
occurs  more  frequently.  This  distortion  is  then  compensated  for  by  weighting  each  event 
by  the  ratio  of  the  probability  that  this  specific  event  would  have  occurred  if  the  true 
probabilities  had  been  used  in  the  simulation  to  the  probability  that  this  same  event  would 
occur  with  the  distorted  probabilities.  Consequently  by  proper  choice  of  the  distorted 
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probabilities  the  number  of  repetitions  can  be  reduced  greatly.  For  instance,  the  mean 
of  a function  Q{x)  is  given  by 

£{Q(JC)}  = J Q(x)dP{x), 

where  P(jc)  is  the  distribution  of  x.  The  mean  of  Q(x)  can  be  estimated  by  selecting  M 
independent  samples  x,  from  P{x)  and  associating  the  probability  1/Af  with  each  event. 
Then  £{Q(x)}  can  be  estimated  by 


(1) 


The  importance-sampling  technique  uses  the  Radon-Nikodyn  derivative  to  express  the 
mean  value  of  Q(x)  by 


£{9(»)|  - / 


where  G(x)  is  a distribution  function.  The  mean  F{Q(x)}  can  be  estimated  by  selecting 
M independent  samples  from  G{x)  and  associating  the  probability  dP(x,)/MdG(x,)  with  each 
event  Q(x,).  Thus  E {Q(x)|  is  estimated  by 


Q(^/) 


dP(Xj) 
dG(Xi)  • 


(2) 


Since  (1)  and  (2)  are  both  unbiased  estimates  of  Q(x),  it  is  possible  to  select  G(x)  so  that 
the  variance  of  (2)  is  less  than  the  variance  of  (1). 


In  our  problem  of  determining  the  threshold  for  a given  Pfa,  when  MTI  samples  are 
noncoherently  integrated,  it  is  necessary  to  estimate  the  distribution  curve 


where 


P(Zj  < T)  = 1 - 


(3) 


(4) 


in  which 


(5) 
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where,  for  a two-pulse  MTI, 


4 = Xij  - Xi.ij  (6) 

and 

y ij  ~ y ij  ~ yi-ij  > c^) 

with  Xjj  and  being  independent  Gaussian  variables  with  zero  mean  and  a variance  of 
a and  P{k)  being  the  noise  power  out  of  a k-pulse  MTI:  P(2)  = 2,  P(3)  = 6,  P(4)  = 20, 
and  P(5)  = 70.  The  straightforward  way  of  estimating  (7)  is  to  generate  Gaussian  samples 
by 

xij  = a(-2  In  sin  27rCy-  (8) 

and 

yjj  = a(-2  In  cos  27ru^j,  (9) 

with  Ujj  and  Vjj  being  independent  random  numbers  uniformly  distributed  on  the  interval 
(0,1).  To  estimate  (3),  M independent  sums  {Zj,j  = l,Af}  are  formed  using  (4)  through 
(7),  and  the  estimated  distribution  is 


M 


P(Z  > T) 


; = 1 


where 


5y  = 1,  Zj  > T, 

= 0,  Zj  < 0. 


Importance  sampling  differs  from  the  previous  procedure  by  generating  samples 

using 


Xij  = a(-2  In  sin  27rCy  (10) 

and 

y,y  = a(-2  In  Uy)'^^  cos  27rCy,  (11) 

where  a > o,  a device  which  yields  more  false  alarms.  Using  (10)  and  (11)  and  (4)  through 
(7),  Af  sums  Zj  are  generated.  Then  the  estimated  distribution  is 


M 


P(Z  > T) 


=-y 
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where 


and 


6y  = 1,  Zj  > T, 

= 0,  Zj  < 0, 


1=2-/' 


27ra2 


e-(xl  * yfj)/2o^ 


+ yf,)l2a^ 


27r«2 


With  use  of  a=  2.0  and  M = 20,000  for  N = 4,  a = 1.7  and  M = 10,000  for  N = 8, 
a = 1.5  and  M = 10,000  for  N = 16,  and  a = 1.3  and  M = 2500  for  N = 32,  threshold  curves 
were  generated  for  two-,  three-,  four-,  and  five-pulse  (binary  weighting)  MTIs  and  are  shown 
in  Fig.  1.  The  reference  curve  for  independent  samples  was  generated  using  detection  curves 
in  Robertson  [3] . 

PROBABILITY  OF  DETECTION 

Since  the  S/N  out  of  the  MTI  is  a function  of  the  target  doppler,  the  doppler 
frequency  where  the  input  and  output  S/N  are  equal  will  bt‘  used.  The  S/N  gain  (or 
loss)  provided  by  the  k-pulse  MTI  is 


1=1 


where  {a^,  i = 1,  ...  , k)  are  the  MTI  coefficients  and  ^ is  the  change  m target  phase 
between  successive  PRFs.  Setting  (12)  equal  to  1 and  solving  for  yields  the  solutions 
02  = 90°,  03  = 103°,  04  = 110.9°,  and  0g  = 116.5°. 

Thus  the  Pd  for  a k-pulse  MTI  and  a given  can  be  found  by  generating  sample 
video  using 

Xjj  = a(-2  In  sin  27ru,y  + A sin  i0^  (13) 

and 

yjj  = a(-2  In  cos  2irVjj  + A cos  10;,,  (14) 
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where  S/N(clB)  = 10  log  (A^I2o^).  By  use  of  (13)  and  (14)  and  (3)  through  (7), 

M = 1024  Zj  values  were  generated  for  each  S/N  and  compared  to  the  appropriate 
threshold.  The  Pp  curves  for  = 10~®  are  shown  in  Fig.  2. 

The  difference  between  the  Pp  curves  for  the  various  MTls  and  the  curve  for  inde- 
pendent pulses  is  the  MTI  noise  integration  loss.  This  loss  is  given  in  Table  1 for  the 
Pd  and  Pfa  values  indicated.  The  loss  appears  to  be  fairly  independent  of  both  N,  the 
number  of  pulses  integrated,  and 


(c)  N = 16  (d)  N = 32 

Fig.  2 — Probability  of  detection  curves  for  N pulses  integrated  with  Pf^  = 10"® 
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Table  1 — MTI  Noise  Integration  Loss  for  Pp  - 0.9 
and  N Noncoherent  Pulses  Integrated 


MTI 

Pulses 

Loss  (dB) 

Average 

Difference 

(dB) 

II 

CO 

II 

N = 16 

II 

CO 

to 

Pfa 

= 10-6 

Two 

1.1 

1.1 

1.1 

0.9 

1.0 

Three 

1.8 

1.9 

1.7 

1.7 

1.8 

Four 

2.2 

2.4 

2.1 

2.1 

2.2 

Five 

2.5 

2.7 

2.3 

2.4 

2.5 

Pfa 

II 

o 

1 

1 

Two 

1 

1.1 

0.9 

0.9 

0.8 

0.9 

Three 

1.8 

1.7 

1.6 

1.5 

1.6 

Four 

2.1 

2.1 

1.9 

1.9 

2.0 

Five 

2.3 

2.5 

2.1 

2.2 

2.3 

COMPARISON  WITH  PREVIOUS  RESULTS 

The  number  of  effective  pulses  integrated  for  a k-pulse  MTI  is  given  [1]  by 


NAk)  = 


N2 


N - 1 

N + 2 ^ (N-j)RlU) 

J = 1 


where  R/jU)  is  the  correlation  coefficient 


Rk(i) 


P(k) 


Thus,  to  find  the  MTI  noise  integration  loss,  the  difference  must  be  found  between  the 
required  S/N  for  and  N independent  pulses.  To  accomplish  this,  a curve  of  S/N  versus 
N for  Pp  = 0.9  and  Pf^  = 10”®  was  generated  using  the  detection  curves  in  Robertson 
[3]  and  is  shown  in  Fig.  3.  From  this  curve  the  MTI  noise  integration  loss  was  calculated 
and  is  shown  in  Table  2.  These  losses  are  about  0.2  dB  higher  than  the  corresponding  losses 
in  Table  1. 
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Fig.  3 — S/N  for  = 0.9  and  = 10"®  as 
a function  of  the  number  of  independent  pulses 
integrated 


Table  2 — MTI  Noise  Integration  Loss  Using  the 
Effective  Number  of  Pulses  Integrated 
for  Pq  =0.9  and  = 10"® 


MTI 

Pulses 

Loss  (dB) 

Average 

Difference 

(dB) 

II 

00 

II 

N = 16 

N = 32 

Two 

1.1 

1.2 

1.2 

1.1 

1.1 

Three 

1.8 

1.7 

1.9 

1.8 

1.8 

Four 

2.4 

2.6 

2.4 

2.3 

2.4 

Five 

2.7 

2.9 

2.9 

2.7 

2.8 

SUMMARY 

MTI  signal  processing  correlates  the  receiver  noise,  and  this  results  in  an  MTI  noise 
integration  loss.  The  losses  for  two-,  three-,  four-,  and  five-pulse  MTIs  are  approximately 
1.0,  1.8,  2.2,  and  2.5  dB  respectively.  The  Pd  for  a given  target  can  be  found  using  the 
following  procedure: 

1.  Calculate  the  input  S/N  (to  the  MTI)  using  the  radar  range  equation; 

2.  Calculate  the  output  S/N  from  the  MTI  using  (12) 

3.  Use  Fig.  2 to  determine  Pp  or  else  assume  all  N pulses  are  independent, 
reduce  S/N  by  the  MTI  noise  integration  loss,  and  find  Pd  from  standard 
detection  curves  such  as  given  in  Robertson  |3). 
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